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Abstract: Reconstructing the hydrological change based on dendrohydrological data has important 
implications for understanding the dynamic distribution and evolution pattern of a given river. The 
widespread, long-living coniferous forests on the Altay Mountains provide a good example for carrying 
out the dendrohydrological studies. In this study, a regional composite tree-ring width chronology 
developed by Larix sibirica Ledeb. and Picea obovata Ledeb. was used to reconstruct a 301-year annual (from 
preceding July to succeeding June) streamflow for the Haba River, which originates in the southern Altay 
Mountains, Xinjiang, China. Results indicated that the reconstructed streamflow series and the 
observations were fitting well, and explained 47.5% of the variation in the observed streamflow of 
1957—2008. Moreover, floods and droughts in 1949—2000 were precisely captured by the streamflow 
reconstruction. Based on the frequencies of the wettest/driest years and decades, we identified the 19% 
century as the century with the largest occurrence of hydrological fluctuations for the last 300 years. After 
applying a 21-year moving average, we found five wet (1724—1758, 1780-1810, 1822-1853, 1931-1967, 
and 1986-2004) and four dry (1759-1779, 1811-1821, 1854-1930, and 1968-1985) periods in the 
streamflow reconstruction. Furthermore, four periods (1770—1796, 1816-1836, 1884-1949, and 
1973-1997) identified by the streamflow series had an obvious increasing trend. The increasing trend of 
streamflow since the 1970s was the biggest in the last 300 years and coincided with the recent 
warming-wetting trend in northwestern China. A significant correlation between streamflow and 
precipitation in the Altay Mountains indicated that the streamflow reconstruction contained not only local, 
but also broad-scale, hydro-climatic signals. The 24-year, 12-year, and 2.2—4.5-year cycles of the 
reconstruction revealed that the streamflow variability of the Haba River may be influenced by solar 
activity and the atmosphere—ocean system. 
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1 Introduction 


Water is a crucial natural resource for providing sufficient food and a sustainable environment for 
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the planet’s human population and for other animals, plants, and microbes (Pimentel et al., 2004). 
Glacial retreat as a result of global warming and water consumption increases as results of 
economic development and population growth is exacerbating water stress in Central Asia (Yao et 
al., 2004; Siegfried et al., 2012). The management and control of water resources based on an 
understanding of the dynamic distribution and evolution pattern of a given river is one of the key 
approaches for addressing water shortages. However, most instrument-based hydrological records 
in Xinjiang, China, span less than 60 years. The short time period of the observed hydrological 
data does not adequately reflect the water resource and river streamflow patterns over a long-term 
period. To address this information gap, we used the tree-ring data to extend our knowledge of the 
hydroclimatic variability back hundreds of years for seasonal or annual timescales. Tree-ring data 
are one of the best sources of proxy streamflow information, because they allow precise dating, 
annual resolution, and comparability with recorded hydro-meteorological data (Pederson et al., 
2001; Woodhouse, 2001; Polacek et al., 2006; Lara et al., 2008; Maxwell et al., 2011; Devineni et 
al., 2013; Shah et al., 2014). 

The Haba River, originating in the Altay Mountains, is one of headstreams feeding the Irtysh 
River, which is the only river flowing to the Arctic Ocean in China. The widespread, long-living 
coniferous forest growing in the Altay Mountains provides a good example for carrying out 
dendrohydrological studies. Li et al. (1989) provided a 290-year length annual streamflow 
reconstruction for the Irtysh River based on tree-ring data. Compared with recent 
dendrohydrological studies (Gou et al., 2007; Yuan et al., 2007; Liu et al., 2010; Yang et al., 2011; 
Bao et al., 2012; Sun et al., 2013), there are deficiencies in existing methods used for streamflow 
reconstruction for the Irtysh River, such as limited sample depth, non-standard measurement, and 
simple stability tests. This study addresses these deficiencies by obtaining more reliable and 
accurate reconstruction using standard dendrochronological procedures. The aims of this study 
are to: (1) construct a tree-ring width based regional composite chronology from larch and spruce 
trees in the mountainous area of the Haba River, (2) reconstruct streamflow variations for the 
Haba River and explore the variations, and (3) assess regional climate signals inhering in the 
developed streamflow reconstruction. 


2 Materials and methods 


2.1 Sample collection 


The Haba River has a length of 214.1 km and origins in the southern Altay Mountains 
(47°52'—49°09'N, 86°06'—87°08'E; Fig. 1), Xinjiang, China. The average annual runoff of the 
Haba River is approximately 2.476x10°? m? and its drainage area covers 7224 km?. The study area 
has a temperate, continental climate. The annual average temperature is 4.7°C, the annual total 
precipitation varies from 100 to 350 mm, and the annual evaporation is 2010.3 mm. 

Two types of long-lived conifer, i.e., larch (Larix sibirica Ledeb.) and spruce (Picea obovata 
Ledeb.) were selected as our materials. The larch trees often grow up to 20 m high with elevations 
ranging from 1200 to 2600 m a.s.l., while these spruce trees often grow up to 35 m high with 
elevations ranging from 1200 to 1800 m a.s.l.. The forest stands are moderately open and the 
canopy densities are low. The soil type at the sampling sites is a Haplic Chernozems, which is thin 
and contains considerable amounts of gravel and rock (FAO et al., 2012). Subalpine shrubs and 
herbs, such as Vaccinium vitis-idaea L., Cotoneaster megalocarpus M. Pop, Lonicera altaica Pall, 
Pyrola rotundifolia L., Aconitum rotundifolium Kar et Kir., Poa supine Schrad., and Carex 
aneurocarpa V. Krecz., are scattered in the understory of the forests (Liu, 1997). 

Information of the three sampling sites, i.e., Qibaakeshi (QBA), Maoyite (MYN) and 
Baolebadai (BLY), is shown in Table 1. To minimize non-climatic effects on tree growth, we 
selected healthy trees with little evidence of fire or human disturbance for sampling. For 
cross-dating purposes, two cores with different directions were extracted from selected trees, 
while only one core was extracted from younger trees. Overall, 161 cores from 84 living trees 
were collected at chest height from the three sites using increment borers (5.15 mm in diameter). 


201712.00325v1 


chinaXiv 


ChinaXiv@ ERAT! 


ZHANG Tongwen et al.: Reconstruction of hydrological changes based on tree-ring data of the Haba... 


86.0°E 86.5°E 87.0°E 87.5°E 


48.5°N a + 48.5°N 


Legend 
Watershed 
— River 
k ET @ Sampling site 
15 «304m E Hydrological station 
E Meteorological station 


a 48.0°N 
86.0°E 86.5°E 87.0°E 87.5°E 


Fig. 1 Map of the Haba River watershed and sampling sites. QBA, Qibaakeshi; MYN, Maoyite; BLY, 
Baolebadai. Abbreviations are the same as in the following tables and figures. 


Table 1 Information of the sampling sites in the Haba River 


, ; Percentage Date of 
Site Latitude Longitude DASA cue Aspect a ers of absent A ee sampling 
B rings (%) E (yyyy-mm-dd) 
405 Larix 
QBA 48.42°N 86.51°E 32/58 1700 S 30 0.328 sibirica 2009-08-01 
(1605-2009) 
Ledeb. 
303 Picea 
MYN 48.46°N 86.70°E 26/52 1150 N 30 0.149 obovata 2010-07-15 
(1708-2010) 
Ledeb. 
305 Picea 
BLY 48.49°N 86.88°E 26/51 1650 NW 25 0.354 obovata 2014-08-20 
(1710-2014) Ledeb 


Note: QBA, Qibaakeshi; MYN, Maoyite; BLY, Baolebadai. 
2.2 Tree-ring width measurements and chronology development 


All the tree-ring cores were naturally dried, mounted, surfaced, and marked following standard 
dendrochronological procedures (Speer, 2010). Each tree-ring width was measured using a 
Velmex Measuring System with a resolution of 0.001 mm (Gou et al., 2007). Cross-dating quality 
control was performed using the COFECHA program (Grissino-Mayer, 2001). The ARSTAN 
program (Cook and Krusic, 2005) was used for detrending and developing tree-ring chronologies. 
The tree-ring width series were detrended with a negative exponential function (a conservative 
function) to remove non-climate signals related to tree age (Liu et al., 2010). For each tree-ring 
width series, the residual chronology was generated by removing the temporal autocorrelation 
using autoregressive moving average time-series models. The pooled model of autoregression 
was reincorporated into a residual chronology to develop an Arstan chronology containing the 
persistent, common, and synchronous signals among for a large proportion of the series from the 
sampled tree cores. Then, three types of chronologies from the sample sites were obtained (Cook 
and Kairiukstis, 1990). Generally, the sample depth within a chronology decreases further back in 
time. This could result in time-dependent variance changes in the earlier part of the chronology 
and a weaker common signal (Woodhouse, 2003). Therefore, the variance in the chronologies was 
stabilized in the chronology compilation process using the method outlined by Osborn et al. 
(1997). Their method uses average correlations of the series in combination with sample size for 
each year to adjust the variance of the changes in the sample size. Subsample signal strength (SSS) 
can be used to assess the adequacy of the replication in the early years of the chronology (Wigley 
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et al., 1984). Therefore, to ensure our study focused only on the period of credible chronology, we 
restricted our analyses to the periods with a SSS of at least 0.85. 


2.3 Hydro-meteorological data 


Monthly mean temperature and precipitation data for 1958—2014 were obtained from the Habahe 
Meteorological Station (48.05°N, 86.40°E; 534 m a.s.l.) near three sampling sites. Climatic data 
were obtained from the China Meteorological Data Service Center (CMDC, 
http://cde.cma.gov.cn/). Monthly streamflow data for 1956-2008 were obtained from the 
Kelatashi Hydrological Gauging Station (48.22°N, 86.42°E; 590 m a.s.l.), which is located in the 
mountain-pass of the Haba River watershed. These hydrological data were provided by the 
Hydrological and Water Resources Bureau of Xinjiang. Figure 2a shows that the highest 
temperatures occur in summer (June—August) and that the highest mean temperature occurs in 
July (22.1°C). A major proportion of the total annual precipitation falls during May-July, 
although there are peaks in both July (21.8 mm) and November (21.4 mm). Much less 
precipitation falls in winter. The climate data recorded in the study area since 1958 show that 
there are significant increasing trends occurring in both temperature (P<0.001) and precipitation 
(P<0.01) (Figs. 2b and c). Figure 2d reveals that streamflow from June to August comprises a 
major proportion of the annual total runoff and peaks in June (204.6 m?/s). The annual streamflow 
fluctuates around the mean value, and the decreasing trend during the period 1956-2008 is not 
significant (Fig. 2e). The correlation coefficients between annual streamflow and annual total 
precipitation (7=0.599, P<0.001, n=51), and between annual streamflow and annual mean 
temperature (r= —0.356, P<0.01, n=51) exceed significance level. 
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Fig. 2 Meteorological data and streamflow data of the Haba River. (a) Monthly precipitation (P) and monthly 
mean temperature (T); (b) Annual precipitation during 1958-2014; (c) Annual mean temperature during 
1958-2014; (d) Monthly streamflow; (e) Annual streamflow during 1956—2008. 


2.4 Statistical analysis 


The developed tree-ring chronologies were decomposed into high- and low-pass components 
using 13 years reciprocal filters to assess the variation patterns in different frequency domains 
(Yuan et al., 2013). The weights of the filters are calculated by the reciprocal filters described by 
Fritts (1976). Pearson’s correlation coefficients were used to quantify the strengths of tree growth 
climate response and climatic and hydrological signals contained in the radial growth of trees. 
After confirming the strongest seasonal relationship between tree-ring width and observed data, 
we used a linear regression model for the reconstruction. Because the length of the 
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instrumental-based records is relatively short for the independent sub-period calibration and 
verification tests, we used the methods of Bootstrap (Young, 1994) and Leave-one-out cross 
validation (LOOCV) (Michaelsen, 1987) to evaluate the statistical reliability of the reconstruction 
model. In the LOOCV analysis, several statistical procedures, including coefficient of efficiency, 
sign test, and product mean test were used to evaluate the coherence between the observed and 
estimated data (Cook et al., 1999). Power spectrum analysis was applied to examine sensible 
periodicities in our reconstruction (Fritts, 1976). Wavelet analysis was used with a Morlet wavelet 
to investigate the periodicity of the reconstructed series and to examine how this periodicity 
changed over time (http://climexp.knmi.nl). Spatial correlation was used to identify coherence 
between our reconstruction and the gridded 0.5°x0.5° Climatic Research Unit (CRU) 
self-calibrating Palmer Drought Severity Index (PDSI) 3.21 dataset (Wells et al., 2004) in a large 
region using the Royal Netherlands Meteorological Institute (KNMI) Climate Explorer. 


3 Results and discussion 


3.1 Statistical characteristics of the tree-ring width chronologies and the development of a 
regional composite chronology 


For further analysis, we used the standard chronologies that contained the common variations 
among the individual tree-core series and retained a high-frequency common variance. This 
variance is presumably in response to climate (Cook, 1985). Because of the low correlations with 
the master dating series, 4 cores (from 3 trees) in the QBA site, 1 core (from 1 tree) in the MYN 
site, and 1 core (from 1 tree) in the BLY site were rejected. Ultimately, 54 cores (from 31 trees) in 
the QBA site, 51 cores (from 26 trees) in the MYN site, and 50 cores (from 25 trees) in the BLY 
site were used to develop ring-width chronologies. Three tree-ring width chronologies with their 
sample depths are shown in Figure 3, and the general statistics of the description of the 
chronologies for a common period analysis from 1900 to 1999 are listed in Table 2. 

The standard deviation (SD) estimates the variability of measurements for the whole series, and 
the mean sensitivity (MS) is an indicator of the relative change in ring-width variance between 
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Fig. 3 Three tree-ring width chronologies (QBA, MYN, and BLY) and one regional ring-width composite 
chronology (WAM) and their sample depths. Solid lines represent the tree-ring width indices and dashed lines 
indicate the sample depths. 
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Table 2 Statistical characteristics of tree-ring chronologies during 1900-1999 


Statistic QBA MYN BLY WAM 
Standard deviation (SD) 0.271 0.299 0.297 0.250 
Mean sensitivity (MS) 0.214 0.240 0.197 0.188 
First-order autocorrelation (AC1) 0.459 0.505 0.636 0.513 
Interseries correlation (trees) 0.360 0.508 0.256 0.245 
Interseries correlation (all series) 0.365 0.515 0.266 0.249 
Mean within-tree correlation 0.541 0.811 0.749 0.658 
Signal-to-noise ratio (SNR) 31.091 49.877 15.621 47.623 
Expressed population signal (EPS) 0.969 0.980 0.940 0.979 
First year SSS>0.85 (core number) 1704 (8) 1773 (5) 1773 (11) 1714 (15) 


Note: SSS, subsample signal strength; WAM, regional ring-width composite chronology. 


consecutive years (Fritts, 1976). Values of SD and MS for the chronologies obtained using the 
regional ring-width composite chronology (WAM) as a result of the combination of tree-ring 
samples are slightly lower than those of the chronologies for a single sample site. The first-order 
autocorrelation (AC1) assesses relationships with previous growth. These values range from 
0.459 to 0.636, which shows that the chronologies contain low-frequency variance generated by 
climate and by tree-physiological lag effects. The interseries correlations of the WAM are also 
relatively lower as a result of the increase in sample depth from the three sampling sites. The 
value of SSS exceeded 0.85, revealing that the threshold corresponds to a minimum sample depth 
of 15 cores for the WAM. Thus, a credible time span for the regional composite chronology is 301 
years (1714-2014). 

After applying the 13-year reciprocal filters, Pearson’s correlation analysis was performed on 
the three sets of data, i.e., the original unfiltered, high-pass filtered, and low-pass filtered data. As 
shown in Table 3, there were significant correlations among the three sampling sites in the 
original, high-, and low-frequency domains over the period 1773-2009. Therefore, we decided to 
combine all the raw tree-ring width data from the three sampling sites (QBA, MYN, and BLY) to 
develop a longer and better-replicated WAM. This was developed using the methods described in 
the above section (tree-ring width measurements and chronology developments). The 
newly-developed regional ring-width composite chronology and its sample depths are shown in 
Figure 3, and the general descriptive statistics were listed in Table 2. 


Table 3 Correlation coefficients among the three sampling sites in original, high-, and low-frequency domains 
during 1773-2009 


Original domain (n=237) High-frequency domain (n=225) Low-frequency domain (n=225) 


QBA MYN BLY QBA MYN BLY QBA MYN BLY 
QBA 1.000 0.431" 0.319* 1.000 0.250* 0.273* 1.000 0.599" 0.363" 
MYN 1.000 0.382" 1.000 0.257° 1.000 0.378" 


BLY 1.000 1.000 1.000 


Note: * indicates the significance at P<0.01 level. 


3.2 Correlation analysis 


A strong biological lag effect was indicated by the high AC1 values for the chronologies from 
single and composite sites ranging from 0.459 to 0.636 (Table 2). Therefore, the monthly 
precipitation (1958-2014), monthly mean temperature (1958—2014), and monthly streamflow data 
(1956-2008) from preceding July to succeeding October were used to assess the effects of 
climatic and hydrological factors on the radial growth of coniferous trees in the study area. As can 
be seen from Figure 4, the results of the correlation indicated that the relationship between 
tree-ring width and precipitation was positive in general, and that a significant correlation 
coefficient was found in preceding July (r=0.411, P<0.01, n=56). However, the response between 
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the radial growth of coniferous trees and temperature was negative, and the WAM was negatively 
correlated with the mean temperature of succeeding June (r= —0.475, n=57), July (r= —0.291, 
n=57), and August (r= —0.262, n=57) at the 95% confidence level. After testing different 
combinations of months, the highest correlation coefficient between the WAM and the 
precipitation was found for the period from preceding April to preceding September (7=0.445, 
P<0.001, n=56); however, the regional composite chronology and mean temperature in 
succeeding June still had the strongest relationship. The monthly streamflow from preceding July 
to succeeding October were also tested for correlations with the regional composite chronology. 
Significant correlations were found from preceding July to succeeding September (except October) 
based on the 95% confidence level, peaking in preceding September, and succeeding February. 
After combination, the highest correlation coefficient between the WAM and July—June 
streamflow was obtained (r=0.690, P<0.001, n=52). Our results show that the influence of 
streamflow on tree-ring growth in the study area is similar to that found in previous 
dendroclimatic studies for North China, such as Picea schrenkiana in the Manas River watershed 
(Yuan et al., 2007), Sabina przewalskii in the Heihe River watershed (Liu et al., 2010; Yang et al., 
2011), Juniperus przewalski in the Yellow River watershed (Gou et al., 2007), and Pinus 
sylvestris var. mongolica in the Yimin River watershed (Bao et al., 2012). 
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Fig. 4 Pearson’s correlations of the WAM and precipitation (a) and mean temperature (b) from the Haba River 
station and the streamflow (c) from the Kelatashi station. Lowercase letters represent months of preceding year, 
and uppercase letters are months of succeeding year. P4-9 indicates precipitation from April to September. To 
indicates mean temperature in June, and S7-6 indicates streamflow from preceding July to succeeding June. 
Dashed lines indicate the 0.05 significance levels. 


3.3 Streamflow reconstruction and stability tests 


We computed the correlation coefficients between the radial growth of larch and spruce trees 
within the study area and various assemblages of months of streamflow for the period 1957—2008 
to determine the appropriate season for reconstruction. To do this, we reconstructed the 
streamflow from preceding July to succeeding June using the regional composite chronology. A 
linear regression model was used to describe the relationship between the tree-ring chronology 
and streamflow. The model was as follows: 

S7-6=31.0814+36.149xWAM, (1) 
where S7-6 is the streamflow from preceding July to succeeding June for the Haba River (m*/s) 
and WAM is the regional tree-ring width composite chronology. The model accounts for 47.5% 
(46.5% after adjusting for loss of degrees of freedom) of the gauged streamflow variance during 
the calibration period 1957—2008 for the WAM. Figure 5a shows that the reconstructed July—June 
streamflow and the observations are well fitting, including the clear low-frequency changes. We 
further compared the first differences (year-to-year changes) of the reconstruction and the 
observations, and obtained a correlation coefficient of 0.595 (P<0.0001, n=51; Fig. 5b). 

The results of the Bootstrap test (100 iterations in the recomputation process) show that all 
values of r, R*, Raq, standard error (SE), F-value (F), P, and Durbin-Watson (DW) are close to 
those of the original regression model 1 (Table 4). Furthermore, the resulting statistics of LOOCV 
are shown in Table 4 (Cook et al., 1999). The results of the low-frequency sign test (S1) and the 
high-frequency sign test (S2), which describe how well the predicted value tracks the direction of 
the actual data, all exceed the 0.01 significant level. The coefficient of efficiency exceeds 0. The 
values of the product mean test (f) are positive, which indicates significant accuracy in the 


JOURNAL OF ARID LAND 


tree-ring estimations (Fritts, 1976). These test results all demonstrate that model 1 is stable and 
reliable and that it could be used for streamflow reconstruction. Consequently, we reconstructed 
the Haba River streamflow from preceding July to succeeding June during the period 1714—2014 
based on model 1, resulting in a mean of 66.74 m?/s and a standard deviation of o=8.89 m?/s (Fig. 


5c). 
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Fig. 5 (a) Comparison between the recorded (blue line) and reconstructed (red line) July—June streamflow for 
the Haba River. (b) Comparison between the first differences (year-to-year changes) of the recorded (blue line) 
and reconstructed (red line) July—June streamflow series. (c) July—June streamflow reconstruction for the Haba 
River since 1714 (solid line). The dashed horizontal lines represent the long-term mean for the period 1714-2014 
and the mean value+lo. Blue and red dots indicate the wettest and driest years, respectively. 


Table 4 Verification results from the Bootstrap and Leave-one-out tests for streamflow reconstruction for the 
Haba River 


Statistic Calibration Bootstrap (100 iterations) Statistic Leave-one-out cross 
mean (ranges) validation (LOOCV) 
r 0.690 0.692 (0.525-0.879) 


Low-frequency 40+/12— (P<0.01) 


nein 0.475 0.484 (0.275-0.773) sign test (S1) 
correlation (R^) 
Adjusted squared multiple 0.465 0.474 (0.258-0.767) 


correlation (R°) High-frequency 37+/14- (P<0.01) 


sign test (S2) 


Standard error (SE) 9.599 9.332 (6.172—12.362) 

F-value (F) 45.318 51.132 (13.535-142.179) Product mean e 
P 0.001 0.001 (0.001-0.001) test (t) 

Durbin-Watson (DW) 1.583 1.222 (0.552-1.953) Coefficient of 0.43 


efficiency (CE) 


3.4 Characteristics of the streamflow reconstruction 


We defined a wet year as >mean+lo (75.63 m?/s) and a dry year as <mean-lo (57.85 m°/s) 
according to the method of division used by Liu et al. (2014). The newly reconstructed July—June 
streamflow series revealed that 51 years could be categorized as ‘wet years’ (accounting for 
16.9% of the total), 47 years categorized as ‘dry years’ (accounting for 15.6% of the total), and 
the remaining 203 years categorized as ‘normal years’ (accounting for 67.5% of the total). The 
values for extreme years and the means of the extreme decades are listed in Table 5. It can be seen 
that the difference between the wettest (1939) and driest (1900) years is 47.36 m°/s, and that the 
difference between the wettest and driest decades (1740s and 1880s, respectively) is 14.99 m?/s. 
Table 5 also shows the long-term means and variation coefficients (Bai et al., 2014) for the 
streamflow reconstruction. The results demonstrate that streamflow in the Haba River was 
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comparatively abundant in the 18" century, and that an obvious decrease in streamflow occurred 
in the 19" century. Thereafter, the streamflow series start to increase and reach a peak in the early 
21% century. The variation coefficient for the 19 century was slightly higher than any other of the 
last 300 years. In addition, the number of anomalies in the 19" century, including the 4 wettest 
and 6 driest years, and the 3 wettest and 5 driest decades, was greater than in any other of the last 
300 years. These findings are consistent with those from a 305-year long early summer 
temperature reconstruction for the southern Altay Mountains (Zhang et al., 2015), which 
demonstrated that there was a greater climatic fluctuation in the 19" century, based on the number 
of extreme years and decades, than in the other two centuries. 


Table 5 Summary of streamflow reconstruction characteristics for the Haba River 


10 wettest years 10 driest years 10 wettest decades 10 driest decades Long-term trend 


Year Mo Year hen Decade Ge Decade an Years nee. Variation coefficient 
1939 94.02 1900 46.66 1740s 74.99 1880s 60.00 1714-1799 67.08 0.131 

1841 93.55 1849 50.13 1840s 72.55 1760s 60.61 1800-1899 66.24 0.136 

1801 90.80 1824 51.18 1990s 71.99 1920s 60.67 1900-1999 66.74 0.134 

1719 89.68 1963 51.18 1800s 71.56 1900s 60.70 2000-2014 68.17 0.132 

1777 88.49 1983 51.54 2000s 71.05 1810s 61.93 1714-2014 66.74 0.133 

1845 86.68 1859 51.58 1830s 70.44 1850s 62.02 

2005 86.21 1880 52.45 1930s 70.39 1870s 63.06 

1751 85.63 1894 52.73 1940s 69.95 1890s 63.87 

1800 84.29 1855 53.02 1750s 69.77 1770s 64.46 


1996 83.35 1768 53.06 1950s 69.39 1970s 64.74 


Unlike the detailed historical documents from southern China, the records regarding 
meteorological disasters in Xinjiang, China, are sporadic and limited as a result of the historically 
low level of human activity. Historical records regarding weather-induced disaster for the Haba 
River watershed are also extremely limited. After the founding of the People’s Republic of China 
(1949), many meteorological stations were established and normal observations were carried out. 
These records allow one of the ten wettest years (1996) and two of the ten driest years (1963 and 
1983) in our hydrologic reconstruction to be compared with the meteorological records for the 
period 1949-2000 (Wen et al., 2006). The results show a good agreement between our 
reconstruction, meteorological records, and meteorological descriptions based on these records 
(Table 6). Furthermore, there are 1776 and 1900 as the extremely wet and dry years respectively 
finding in the streamflow reconstruction for the upper Irtysh River (Chen et al., 2016) and 
precipitation reconstruction for the southern Altay Mountains (Chen et al., 2014). A very severe 
drought period during the 1920s occurring in a wide area of northern China is mentioned in 
historical records (Shi et al., 1991) and tree-ring based hydro-climatic reconstructions (Liang et 
al., 2006; Gou et al., 2007). This consistent drought is also recorded as the third driest decades in 
the Haba River streamflow series (Table 5). Thus, our reconstructed hydrological series can 
capture the signals of flood or drought disasters for the Haba River watershed. 


3.5 Linkages to the recent warming-wetting trend 


The decadal variability was highlighted by applying a 2l-year moving average to the 
reconstruction, as shown in Figure 6, from which five wet and four dry periods can be 
distinguished. The wet periods (above the mean value of the reconstruction) were 1724-1758 
(average value: 69.03 m?/s), 1780-1810 (68.09 m?/s), 1822-1853 (69.50 m?/s), 1931—1967 (69.11 
m/s), and 1986-2004 (70.24 m/s). The dry periods (below the mean) were 1759-1779 
(63.67m°/s), 1811—1821 (65.16 m?/s), 1854-1930 (63.49 m/s), and 1968-1985 (64.65 m°/s). 
Figure 6 reveals that there were 4 periods, i.e., 1770-1796 (lasting for 27 years), 1816—1836 (21 
years), 1884-1949 (66 years), and 1973-1997 (25 years) with an increasing trend existing in the 
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streamflow reconstruction in the low-frequency domain. 


Table 6 Wettest and driest years of in the hydrological reconstruction of the Haba River watershed in 
comparison with meteorological records for the period 1949-2000 


Type Year Short description of flood or drought disaster 


1. Some floods occurred in the Altay region in July and August, 1995. 
2. Some floods occurred in the Altay region in June, July, August, and September, 1996. 
1. Great drought occurred in the whole Xinjiang in 1962. There was 50% streamflow of the Altay 
River, and the streamflow of the Burqin River reduced 40%. People and livestock lacked adequate 
Driest year 1963 drinking water. 
2. Drought occurred in the Altay, Tacheng, and Ili regions in 1963. Less rain and snow in spring 
led to drought in the Altay region. 
1. Drought occurred in the Altay, Tacheng, Ili, and Karamay regions in 1982. There was 10% 
production reduction in the irrigated land, while there was no harvest in dry land in the Burgin 
region. 

1983 2. Infrequent drought occurred in the whole Xinjiang in 1983. Extraordinary drought occurred 
due to the rare snow in winter and less rain in the following season in the Burqin region. The 
streamflow of the Burgin River was 13 m°/s in April, and the streamflow of the Irtysh River was 
only 0.5 m3/s on 3 May. 


Wettest year 1996 
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Fig. 6 Variations in the streamflow reconstruction for the low-frequency domain. The black line shows data 
smoothed with a 21-year low-pass filter to reflect the long-term fluctuations. The red dashed lines represent trend 
lines for the periods with increasing streamflow. The horizontal dashed line represents the long-term mean for the 
period 1714-2014. The blue bars indicate the wet periods; the red bars indicate the dry periods. 


Generally, more precipitation leads to more streamflow as a result of increased runoff, while 
higher temperatures lead to less streamflow because of stronger evaporation. However, rivers in 
the arid and semi-arid areas may have large streamflow despite increasing temperature (Chen and 
Wang, 2004). Based on daily meteorological data, Li and Jiang (2007) found that an abrupt 
change in precipitation occurred (from less to more) over northern Xinjiang in 1984. A warming 
trend since the 1980s has frequently been observed in records obtained from instrument-based 
monitoring (Li et al., 2006). Shi et al. (2007) also indicated that since 1987, the climate has 
changed from the warm-dry to the warm-wet in Xinjiang, China. As a result of the higher 
coefficient of X and R? in the tendency equation, the increasing trend from 1973 to 1997 in the 
newly reconstructed streamflow series for the Haba River changed more dramatically than at any 
other time in the last 300 years. The results confirmed that the streamflow of the Haba River, and 
other small drainage basins in the Xinjiang, had an obvious increasing trend in response to the 
recent warming-wetting trend in northwestern China. 


3.6 Regional/large-scale climate signals inhering in the streamflow reconstruction 


Dai et al. (2004) revealed that basin-averaged annual PDSI covaried closely with streamflow for 
seven of the world’s largest rivers and several smaller rivers. Thus, PDSI can be regarded as a 
good proxy of both surface moisture conditions and streamflow. A significant correlation (7=0.583, 
P<0.01) was found between a 25l-year long drought reconstruction for the northern Altay 
Mountains and the annual streamflow for the Irtysh River (Chen et al., 2015). The spatial 
correlation was applied to evaluate the regional significance of tree-ring width based streamflow 
reconstruction. The streamflow reconstruction correlates (coefficient>0.5) with the gridded 
July-June PDSI data for 1956-2012 in a large area between approximately 43°-51°N and 
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79°—90°E, with the highest correlations (coefficient>0.6) occurring in the northern Junggar Basin 
and the entire Altay Mountains (Fig. 7). 
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Fig. 7 Spatial correlation of the streamflow reconstruction for the Haba River with the gridded Palmer Drought 
Severity Index dataset 


The correlations between the tree-ring width chronologies in the study area and precipitation 
were not high. However, the streamflow reconstructions in the previous studies are limited. 
Therefore, a July—June precipitation reconstruction for the Altay Mountains (an explained 
variance of 65.5%; Chen et al., 2014) was selected to examine whether there is a connection 
between our reconstruction and other reconstructions. The correlation coefficient for the 
streamflow reconstruction and the precipitation reconstruction is 0.54 (P<0.0001, n=185). After 
employing a 13-year reciprocal filter, both the tree-ring width based reconstructions were 
decomposed into high- and low-pass components for comparison. The correlation coefficients in 
the high- and low-frequency domains are 0.44 and 0.59, respectively. Figure 8a shows that the 
streamflow reconstruction and the precipitation reconstruction corresponded well during the 
period 1825-2009. Furthermore, most of the high and low values in the two reconstructions 
appeared in similar years in the high-frequency domain (Fig. 8b). In addition, the variation modes 
and long-term trends of two reconstructions were roughly synchronous in the low-frequency 
domain (Fig. 8c). These results confirmed that broad-scale regional climatic elements could affect 
the streamflow variations of the Haba River. 

The results of the power spectral analysis over the entire reconstruction period (1714-2014) 
revealed that significant periodicity is further indicated at a frequency of 4.5-year (90%), 3.2-year 
(95%), and 2.2-year (90%; Fig. 9a). The temporal characteristics of the different cycles were also 
evaluated using wavelet analysis and are illustrated in Figure 9b. The wavelet analysis indicates a 
robust approximate 24-year cycle from 1880 to 1950. In contrast, a significant (although less 
robust) power with an approximate 12-year cycle is detected from 1780 to 1810 and from 1940 to 
1980. The 12-year cycle in our reconstructed streamflow series suggested a possible important 
connection with the 11-year Schwabe cycle of solar activity (Nagovitsyn, 1997). The 24-year 
cycle roughly corresponded to the 22-year Hale cycle—the double harmonics of the 11-year cycle. 
The inter-annual cycles (4.5-year and 3.2-year) in our reconstruction fall within the range of the 
El Nifio-Southern Oscillation (ENSO) (Allan et al., 1996). Luo (2005) also indicated that there is 
an approximate 3-year cycle appearing in the precipitation variability for the northwestern China. 
Furthermore, the above cycles contained in the streamflow reconstruction for the Haba River can 
often be found in earlier dendroclimatic studies of the Altay Mountains. The 11-year and 2-year 
cycles contained in the 310-year early summer temperature series for northeastern Kazakhstan 
have been examined using power spectral analysis by Shang et al. (2011). Zhang et al. (2015) 
identified the frequencies of 1l-year and 2.0—2.7-year cycles in a summer precipitation 
reconstruction of the southern slope of the Altay Mountains recorded using tree-ring 5'°C. Jiang et 
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al. (2016) also revealed that there were 24.3-year, 3.2-year, and 2.1-year cycles in the July—June 
precipitation variation for the upper Irtysh River Basin in the last 291 years. 
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Fig. 8 Comparison of the streamflow reconstruction for the Haba River and the precipitation reconstruction for 
the Altay Mountains using 13-year reciprocal filters in the (a) original, (b) high-, and (c) low-frequency domains 
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Fig. 9 Power spectrum (a) and wavelet power spectrum (b) of the reconstructed streamflow series for the Haba 
River. The blue line represents the 90% confidence level, and the red line represents the 95% confidence level. 
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4 Conclusions 


A 301-year long regional composite tree-ring width chronology was developed from 32 larch and 
52 spruce trees sampled from three sites in the Haba River watershed. Because there was a high 
correlation coefficient between instrument-based hydrological data and the radial growth of 
sampling trees, it was possible to develop an annual streamflow reconstruction (from preceding 
July to succeeding June) for the Haba River using the regional composite chronology. The 
reconstruction matched the observed data well, and precisely captured flood and drought events 
occurring in 1776, 1900, 1963, 1983, and 1996 mentioned in the meteorological records and other 
tree-ring based hydro-climatic reconstructions. The 19" century had the 4 wettest and 6 driest 
years, and the 3 wettest and 5 driest decades, and had the largest hydrological fluctuations for the 
target river in the last three centuries. After applying a 21-year moving average, five wet periods 
occurred during the periods 1724-1758, 1780-1810, 1822-1853, 1931-1967, and 1986-2004, 
while four periods 1759-1779, 1811—1821, 1854-1930, and 1968-1985 were relatively dry in the 
low-frequency domain. Furthermore, four obvious increasing trend periods (1770-1796, 
1816-1836, 1884-1949, and 1973-1997) were found in the reconstruction. The period since the 
1970s has changed more dramatically than any other period in the last 300 years as a response to 
the recent warming-wetting trend for the northwestern China. The spatial correlations between the 
streamflow reconstruction and the gridded PDSI dataset revealed a significant positive correlation 
in the area of the northern Junggar Basin and the entire Altay Mountains. The significant 
correlation between our streamflow reconstruction and a precipitation reconstruction for the Altay 
Mountains in the original, high-, and low-frequency domains indicated that the streamflow 
reconstruction contained not only local but also broad-scale hydro-climatic signals. The 24-, 12-, 
and 2.2-4.5-year cycles in the reconstruction suggested that the streamflow variability of the 
Haba River may be influenced by solar activity and the atmosphere—ocean system. 
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